Numeric Experiments in Relativistic Thermodynamics: 
A Moving System Appears Cooler 
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Abstract 



In this paper we simulate a two dimensional relativistic ideal gas by implementing a relativistic elastic 
binary collision algorithm. We show that the relativistic gas faithfully obeys Jiittner's speed distribu- 
tion function. Furthermore, using this numeric simulation in conjunction with the energy equipartition 
theorem for a relativistic gas, we conclude that a moving system appears cooler. 

cr: 

;-i ! 1 Introduction 

There is a long standing controversy on the rendering of the fundamental thermodynamic concepts in the 

fT^ I special-relativistic form. 

^ I In 1907 Mosengeil, Planck, and Einstein [Hill [3], independently showed that if a thermodynamic system 

^^ . with temperature T is at rest in the inertial frame K, then its temperature T' measured by an inertial frame 

iJ I K' moving with speed u relative to K, is T' — r/7 where 7 = (1 — u^)~^/^ (in c = 1 system of units). 

^.f-\ . Therefore, a moving system appears cooler. However, a completely opposite result, T' — Tj, was proposed 

later by Blanusa [5 , Einstein (in a little known letter to von Laue) [13], Ott [7], Arzelies [TU] and many 
others. Then, Einstein [13 (in another letter to von Laue), Landsberg [8], and others, pondered whether 

QQ ' the temperature should be a relativistic invariant. Few years later, Balescu [5] explored the possibility that 

(^ . all relativistic transformations are in fact equivalent via some sort of "gauge" transformations. This list is 

manifestly incomplete. The relativistic thermodynamics, in its special-relativistic formulation, has a long 
and convoluted history. A detailed discussion about this curious episode of modern physics is beyond the 

K> I scope of this paper. The interested reader should consult for example [TOl [12l [13l [M] and the references 

therein. 

It was often argued that because of the high temperatures, and the relativistic speeds involved, there 
are no experimental results that could settle the controversy on the relativistic thermodynamics. However, 
by deploying a realistic computer model, we can actually "experimentally" decide on this question. Using 
a numeric model of a two dimensional relativistic ideal gas, in this paper we show that a) the relativistic 
gas faithfully obeys Jiittner's speed distribution [4], and b) defining the temperature via the relativistic 
equipartition theorem, we show that the numeric simulations strongly favor the Mosengeil, Planck, and 
Einstein relativistic formula. Thus, we have a compelling numerical indication that moving objects should 
appear cooler. 
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2 Computer Simulations 

We consider a two-dimensional ideal relativistic gas whose particles experience elastic collisions. The gas is 
enclosed inside a rectangular container, with perfectly reflecting walls, which conserve the energy-momentum 
of the colliding particles. Throughout this paper we will use a system of units where the speed of light, the 
mass of the particles, and the Boltzmann's constant are all taken equal to one c = l,m = l,kB = 1. 

The thermodynamic system is at rest in the lab when the container is at rest, and the average velocity 
of its particles is zero {v ) = 0. If the gas is in thermal equilibrium at temperature T — 1//3, its speed 
distribution is given by Jiittner's formula [H [6] 

f{v)=v^texp{-f]-f,)/Z , (1) 

where 7t, = (1 — v^)^^'^ is the Lorentz factor, and Z = e~^{l + P)/ 0^ is a normalization constant, chosen 
such that /p f(v)dv — 1. Then, the average energy is given by {E) = /p Ef{v)dv = (/3^ + 2f3 + 2)/(/3^ -I- /3) . 
We will use this information to find the reciprocal of the temperature of the system at rest, from measuring 
its average energy. We get 

P={2-{E) + ^{Ey+4{E)-4)/i2{E)-2) . (2) 

2.1 Numeric results 

For the experimental setup we considered 100, 000 identical particles, which initially, were randomly dis- 
tributed inside a rectangular box, and given the same speed in arbitrary directions. During one time step 
the algorithm checks for possible collisions with a) the wall of the container (the box has the edges parallel 
with X and y axes), and b) with another particle. If none occurs, the particle advances one time step with 
the same velocity. 

a) Collisions particle-wall: If the particle hits the left or right wall, it reflects around the a;-axis. If it hits 
the top or the bottom wall, the particle reflects around the y-axis. 

b) Collisions particle-particle: Point-like particles have practically no chance to collide if only contact, zero- 
range interactions are present. To bypass this difficulty, we define an "interaction area" around each particle 
as follows. If {xi,yi) are the coordinates of particle i as measured by an observer at rest with respect to 
the container, and if particle j arrives at (xj,yj) such that {xi — XjY + (j/i — 2/j)^ < R^i then the two 
particles "collide". The distance R is chosen to be much smaller than the linear dimensions of the container. 
The algorithm "locks" the two particles when colliding, thus ignoring the possibility of a third particle to 
participate in collision. This simplifies the simulation and does not introduce errors greater than a one 
percent. Knowing the "in" states of the particles entering the collision we determine the "out" states by 
using only the conservation of the energy-momentum. However, in this step, which is solved in the center of 
mass of the two particles, one has to choose an "interaction angle" 9 for the outgoing particles, as illustrated 
in figure ([1]). Since particles are modeled as zero-size points experiencing only elastic collisions, 9 is not 
meaningfully defined. Therefore, we use a pseudo-random number generator to assign a random value to the 
"interaction" angle 6 . 

We monitored the average energy per particle, and considered that the system reached thermal equilib- 
rium when the average energy per particle fluctuated by less than 0.001%. From the onset of the experiment, 
the equilibrium was usually achieved after 40, 000 time steps. After the system reached thermal equilibrium, 
we determined its temperature by means of eq.®. By plugging the measured value back into eq.((T]) we 
were able to compare the theoretical curves with the histograms of the speed distribution obtained from 
our computer modeling. In figure ([2]) we illustrate the experimental versus the theoretical results for the 
2-d relativistic gas. The figure collects data from three different experiments: low temperature (with the 
corresponding reciprocal temperature /3 = mc^ /ksT = 49.46), intermediate temperature (/3 = 7.25), and 
respectively high temperature (with the corresponding reciprocal temperature /? = 2.00). The agreement 
with Jiittner's distribution is remarkable in all three cases. 




Figure 1: The interaction as seen in the center of mass of the two particles. Here gi and (73 are the momenta of 
the incoming and respectively outgoing particle 1. Similarly, 92 and (74 are the momenta of the incoming and 
respectively outgoing particle 2. The angle 9 between incoming and outgoing particles is chosen randomly. 
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Figure 2: Jiittner's speed distribution curves vs. the experimental speed distribution for a relativistic 2-d gas 
at reciprocal temperatures f3 — 49.46, /3 = 7.25, and f3 = 2.00 respectively. The area under each distribution 
curve is equal to the total number of particles, N = 100, 000. 

2.2 The relativistic energy equipartition 

Using the Jiittner distribution function llj we can prove that for the system at rest, we obtain [12l ?] 

1//3 = T = {p^ V^) = {Py Vy) . (3) 

We checked experimentally the validity of formula ^ for all three reciprocal temperatures (3 = 49.46 , 7.25 
and 2.00. The results are summarized in table 1, where we have denoted fix = (PxVx) and (3y = (pyVy). 
Note that the numeric simulations show a good agreement with the theory. 



2.3 The relativistic transformation of temperature 

The system is in translation with the speed u with respect to the lab if the average velocity of all its particles 
measured by a stationary observer in the lab is (v) — u. Without loss of generality we will consider the 
system moving along the x-axis. Extending the equipartition theorem to the moving system K'^ we obtain 



i/p'^T' = {p:,iv^-u)) = {p'v:) 



(4) 



Pv 



2.0000 


2.0120 


2.0147 


7.2476 


7.27541 


7.2487 


49.4561 


49.6334 


49.3251 



Table 1: Energy equipartition for the relativistic gas 



For a complete proof of eq. @ the reader should consult the reference p[2 • Equation ([U is a good tool to 
decide which relativistic transformation temperature formula is consistent with this experiment. As noted 
in [llj . because the thermodynamic system is an extended object, the output of the measurements depends 
whether the measurement is performed at constant time in the rest frame or at constant time in the moving 
system. The extended geometrical dimensions of the thermodynamic system implies that a constant time 
in the rest frame corresponds to a set of times t' in the moving system and vice-versa. We measured the 
particles in the container at constant time t with respect to its rest frame. 



Planck's formalism gives /3p;^„^^ 
is the usual Lorentz factor. 



— Pju, while Ott's formalism gives Pq^^ = /3/7„, where ^^ — (1"""^) ^^^ 
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(a) The relativistic transformation of temperature for 
a relative low rest reciprocal temperature /3 = 49.46. 
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(b) The relativistic transformation of the temperature 
for relative medium (/? = 7.25) and high (/3 = 2.00) rest 
reciprocal temperatures. 



Figure 3: Note that for the boost speeds u = 0.2, 0.4, 0.6, and 0.8 the agreement between (3J^,(3y and Ppianck 
is remarkable. We observe that with the increase of the temperature, P^ shifts slightly below the values of 
l3y. In contrast, /Jq^j is clearly diverging from the experimental trend. 



Choosing the same rest reciprocal temperatures as before (/3 = 49.46, 7.25 and 2.00), we tested eq. ^ 
for the boost speeds u = 0.2, 0.4, 0.6, and respectively 0.8. The results are presented in figures (3(a) I and 
( |3(b)[ ). The robustness of the numeric algorithm was checked by repeating the experiments with different 
form factors for the rectangular box. 



3 Summary 

In this paper we solved the non-trivial problem of relativistic two-dimensional molecular dynamics. We 
showed that a relativistic ideal gas faithfully satisfies Jiittner's speed distribution function. Using as guidance 
the relativistic equipartition theorem, we showed experimentally that the numeric simulations strongly favor 
the Mosengeil, Planck, and Einstein relativistic formula. Thus, moving objects appear cooler. 



Acknowledgements 

We thank Asim Gangopadhyaya, Ovidiu Lipan and Mircea Pigli for valuable comments and discussions. 

References 



[1] 

[2] 

[3] 

[4] 

[5] 
[6] 

[7 
[8 
[9 

[lo; 

[11 

[12 
[13 
[14 



K. von Mosengeil, Ann. Physik, Vol. 327, Issue 5, 867 (1907). 

M. Planck, Sitzber. Preuss. Akad. Wiss. Berlin 542 (1907). 

A. Einstein, Jahrb. Rad. E. 4, 411 (1907). 

F. Jiittner, Ann. der Physik, Vol. 339, Issue 5, 856 (1911). 

D. Blanusa, Glasnik Mat.DFiz. i Astr., Ser. II 2, 249 (1947). 

J. L. Synge, The Relativistic Gas, Interscience Publishers, Inc., New York, (1957). 

H. Ott, Zeit. ffir Physik 175 70 (1963). 

P.T. Landsberg, Proc. Phys. Soc.89 1007 (1966). 

R. Balescu, Physica 40 309 (1968). 

H. Arzelies, Thermodynmique Relativiste et Quantique, Gauthier-Villars Eds., Paris (1969). 

C. K. Yuen, Am. J. Phys. 38: 246-252 (1970) 

A. Guessous, Thermodynamique Relativiste, Gauthier-Villars Eds., Paris (1970). 

C. Liu, Einstein and Relativistic Thermodynamics, Brit. J. Hist. Sci. 25 185 (1992). 

P.T. Landsberg and G.E.A. Matsas, Physica A 340 92 (2004). 



